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Abstract 

By using the grand canonical Monte Carlo simulations on spherical surfaces with 
two fixed vertices separated by the distance L, we find that the second-order phase 
transition changes to the first-order one when L is sufficiently large. We find that 
string tension a ^ in the smooth phase while a —* in the wrinkled phase. 
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1 Introduction 



In the past two decades, a considerable progress has been made [1,2,3,4,5] 
on understanding the phase structure of Helfrich [6] and Polyakov-Kleinert 
[7,8] model, which is an elastic surface model with bending rigidity denoted 
by b. The surface fluctuation can be viewed as a second order phase transi- 
tion between the smooth phase at b —>■ oo and the crumpled phase at b — > 
[9,10,11,12,13,14]. Numerical studies have also been concentrated on the phase 
transition [15,16,17,18,19,20,21,22,23,24,25,26]. 

The scaling of the string tension a was first investigated by Ambjorn et. al 
[27,28]. It has been recognized that a vanishes at the critical point of the 
phase transition [4]. The canonical Monte Carlo was used to extract informa- 
tion about a. However, a is originally defined in the grand canonical ensemble. 
Therefore, it is interesting to use grand canonical Monte Carlo simulations in 
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order to see the scaling of a more convincingly, although it has been demon- 
strated that the canonical MC can simulate the scaling of a [28]. 

It was recently reported that the phase structure of the model depends on 
the choice of the integration measure Y\i<l?dXi, where is the co-ordination 
number of the vertex i [25]. From a conformal field theoretical viewpoint, this 
a is believed to be la = 3 [29,30,31,32]. On the other hand, is considered as 
a volume weight of the vertex i in the integration dXi. Hence, it is possible to 
extend a to continuous numbers by assuming that the weight can be chosen 
arbitrarily. Therefore, it is interesting to see the dependence of a on the phase 
transitions which can be controlled by the parameter a. 

In this article, we would like to investigate the scaling property of a. It will be 
confirmed by using the grand canonical MC that the string tension vanishes 
at the continuous transition. Moreover, we would like to show a numerical 
evidence of the nonvanishing string tension in the smooth phase, which is 
distinguished from the wrinkled phase by the discontinuous transition at finite 
negative value of a. 

We speculatively comment on why the result of nonvanishing string tension 
could be a relevant one. It is possible to consider that the nonvanishing string 
tension is connected to two interesting problems. The first is the problem of 
quark confinement, which is rather a problem in mathematical physics. The 
linear potential V(L) ~ L assumed between quark and anti-quark separated 
by the distance L gives the finite string tension, which is compatible with our 
result of nonvanishing string tension. The second is the conversion of external 
forces into an internal energy and vice versa in real physical membranes, and 
is a rather practical problem. If the model in this Letter represents properties 
in some real membranes, our result implies a possibility of such conversion. 



2 The model 

A sphere in R 3 is discretized with piecewise linear triangles. Every vertex is 
connected to its neighboring vertices by bonds, which are the edges of triangles. 
Two vertices are fixed as the boundary points separated by the distance L. 

The Gaussian energy Si and the bending energy S 2 are defined by 

s ± = x; (x, - x 3 f , s 2 = J2 (i - cos e<) , (i) 

(ij) » 

where J2(ij) is the sum over all bonds (ij), and Q{ in S 2 is the angle between 
two triangles sharing the edge i. 
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The partition function is defined by 



Z(b, fj,, a; L) = ]T ]T / J] dXi exp [S(X, T)} , (2) 

N T i=l 

S(X, T, N) = S 1 + bS 2 - fiN - a ]T log ft , 

where J2t denotes the sum over all possible triangulations T, N is the total 
number of vertex. It should be noted that the chemical potential term — fiN 
and the co-ordination dependent term —a J2i log (ft are included in the Hamil- 
tonian. The expression S(X, T, N) shows that S explicitly depends on the 
variables X, T and N. The coefficient b is the bending rigidity, and ji is the 
chemical potential. Z depends on b, /i, a, and L. The surfaces are allowed to 
self-intersect and hence phantom. 

We expect 

Z(b, /i, a; L) ~ exp(— aL) (3) 



in the limit L — > oo [28]. Then, by using the scale invariance of the partition 
function, we have [4,28] 

a= m^m, ( 4) 



where (^i) and (N) are the mean values of Si and N. 

The specific heat, which is the fluctuation of S 2 , is defined by 
C S2 = {b 2 / (N))(d 2 log Z/db 2 ), and is calculated by using 

Cs 2 = ^((S 2 -(S 2 )) 2 ). (5) 



3 Monte Carlo technique 



X is updated so that X' = X + 5X, where the small change 5X is made at 
random in a small sphere centered at X. The radius 5r of the small sphere is 
chosen to maintain the rate of acceptance r x for the X-update as 0.5 < r x < 
0.55. 5r is defined by using a constant number e as an input parameter so 
that 5r = e (I), where (/) is the mean value of bond length computed at every 
250 MCS (Monte Carlo sweeps). It should be noted that Sr is almost fixed 
because (I) is constant and unchanged in the equilibrium configurations. 
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T is updated by nipping a bond shared by two triangles. The bonds are la- 
beled by sequential numbers and chosen randomly to be flipped. The rate of 
acceptance r T for the bond flip is uncontrollable, and the value of r T is about 
30% < r T < 40%. iV-trials for the updates of X and iV-trials for T are done 
consecutively and these make one MCS. 

N is updated by both adsorption and desorption. In the desorption, a vertex 
is chosen at random, and then a bond that is connected to the vertex is 
chosen at random so that the two vertices at the ends of the bond unite 
and become a new vertex. In the adsorption, a triangle is chosen at random 
just like a bond chosen in the desorption, and a new vertex is added to the 
center of the triangle. As a consequence, the Euler number (=2) of the surface 
remains unchanged by the adsorption/desorption. The acceptance rate is 
uncontrollable as well as r T , and the value of r N is about 55% < r N < 65% in 
our MC. 

In the adsorption of a vertex, the corresponding change of the total en- 
ergy AS = 5'(new) — S'(old) is calculated. The adsorption is then accepted 
with the probability Min[l, exp (— AS) /(N + 1)]. In the desorption, AS = 
«S(new) — S'(old) is calculated by assuming that one vertex is removed. The 
desorption is then accepted with the probability Min [1, iVexp (— Zi,!?)]. The 
adsorption/desorption are tried alternately at every 5 MCS. 

We use surfaces of size iV~500, iV~i000, and iV~i500. N depends on both 
ji and a which is fixed to a = —5.5, and iV is almost independent of L. There 
is no a piori reason for choosing a = —5.5, although we have confirmed that 
the phase transition occurs at a~— 5.5 in the model with the fixed center of 
surface. The values of \i are chosen so that iV~500, iV~1000, and iV~1500. 
The diameter L (N) of the spheres at the start is fixed so that J2i^i— 3iV/2, 
where k is the length of the bond i. As a consequence, L (N) becomes 

L (N) oc y/N. (6) 



We use two kinds of L for each L (N) so that L—1.5L (N) and L = 3L (N). 
The distance of the boundary points is increased from L (N) to L in the first 
5xl0 6 MCS. 

It should be noted that both L = 1.5L (N) and L = 3L (N) becomes oo in 
the thermodynamic limit iV — > oo because of Eq.(6). Therefore, a defined by 
Eq. (3) can be extracted from these values of L at sufficiently large N. 
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Fig. 1. C S2 vs b at (a) a = -5.5, L = l.5L (N) and (b) a = -5.5, L = 3L (N). 
/i = 15.9 (0)> M = 15.5 (A), and ^ = 14.8 (□) respectively correspond to iV~1500, 
iV~1000, and iV~500. 

4 Results 



The specific heat Cs 2 at a = —5.5 is plotted against b in Figs. 1(a) and 1(b). 
The number of vertices is N ~ 500(D), iV ~ 1000(A), and ~ 1500(O) in 
each figure. We find that the peak value C!? ax grows with increasing N. An 
interesting point to emphasize is that the phase transition is strengthened 
when L is increased from L = 1.5L (N) to L = 3L (N). 
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(a) a=-5.5, L=1.5L 
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(b) a=-5.5,L=3L 
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Fig. 2. C!? ax vs (N) in log-log scale at (a) a = -5.5, L = 1.5L , and at (b) 
a = — 5.5,L = 3Lo- 

To see the order of the transition, we plot C<? ax vs (N) in Figs. 2(a) and 2(b) in 
log-log scale. We clearly find in both of the figures that CJ? ax scales according 
to <^ 2 ax ~ N", and we have 



v = 0.265 ± 0.025 [a = -5.5, L = l.5L (N)} , 
v = 0.822 ±0.182 [a — —5.5, L — 3L (N)] . 



(7) 
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The value v — 0.265(25) suggests that the phase transition is of second order at 
a— — 5.5 when L = 1.5L (N). On the other hand at a = — 5.5, L = 3L (N), the 
value z/ = 0. 822(182) is still slightly smaller than v — \. Nevertheless, the result 
implies that the model undergoes the discontinuous transition at a = —5.5 
when L = 3L (N). 
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Fig. 3. Variation of (S2)/(Nb) against MCS at (a) b= 1.86, L = 1.5Lq (at the critical 
point of the continuous transition) and at (b) 6= 1.86, L = 3Lq (at the discontinuous 
transition point), where a = —5.5, fi = 15.9. 

To see the difference between the transitions shown in Figs. 2(a) and 2(b), we 
depict (S2)/{Nb) against MCS in Figs. 3(a) and 3(b). These were obtained 
at b = 1.86, L = 1.5L ,/i = 15.9 and at b = 1.86, L = 3L , \i = 15.9, where 
N B is the total number of bond. It is possible to see in Fig. 3(b) that there 
are two distinct states which differ in values of S2', one is characterized by 
(S 2 )/(N B )~0.2 and the other by (S 2 )/(N B ) -0.25. 

Large number of MCS was needed to obtain v in Eq. (7). 9.5 x 10 s MCS 
were done for iV ~ 1500 surfaces at the vicinity of the transition point when 
L = 3L (N). Relatively small number of MCS was done when L = 1.5L (A^). 

The dependence of a on b is plotted in Figs. 4(a), 4(b), where the symbols 
(O, A, □) in the figures correspond to those in Figs. 1(a), 1(b). We immedi- 
ately understand from the figures that the dependence of a on b in Fig. 4(a) is 
completely different from that in Fig. 4(b). a decreases with increasing b when 
L = 1.5L (N). On the contrary, a increases with increasing b when L = 3L (N). 
Moreover, we find in Fig. 4(b) that o"^o" (const) in the smooth phase close at 
6 c (/i, a), and that <jq is independent of N, where b c (fj,, a) denotes the value of 
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Fig. 4. a ys b obtained at (a) a = —5.5, L = 1.5Lq and at (b) a = — 5.5, L = 3Lq. 

b where Cs 2 has the peak as shown in Fig. 1(b), and denotes that b c depends 
on fi and a. Figure 4(b) also shows that a rapidly changes against b at 6 c (/i, a) 
when N (or jj) increases. 

In order to show the dependence of a on (N) more clearly, we introduce the 
reduced bending rigidity 



A = 



b c (n,a) 



1. 



(8) 



Then, the transition point b c (fi,a) is represented by A = 0, the smooth phase 
at b>b c (/i,a) by A>0, and the wrinkled phase at b<b c (/i,a) by A<0. 
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Fig. 5. a vs (N) in log-log scale obtained at (a) a = —5.5, L = 1.5Lq and at (b) 
a = — 5.5, L = 3Lq. 

Figures 5(a), 5(b) show log-log plots of a against (N) obtained at A = —0.1(A), 
A = 0(O)> and A = 0.1(D). The straight lines in each figure denote the scaling 
property of a such as 



a oc N~ K 



(k > 0). 



(9) 
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It should be emphasized that the scaling in Eq. (9) is compatible with a oc 
(L/N) s in [28], since L oc L (N) oc y/N as described in Eq. (6). We find in 
Figs. 5(a) that er ^0(N — >oo), which is the scaling property at the continuous 
transition in [28]. 

On the contrary, we clearly see in Fig. 5(b) that 

a — > 0" O ( const), [b > b c (fi, a)] . (10) 



This corresponds k — > in Eq. (9), and is the main result of this Letter. 
Moreover, we see a—>0(N — >oo) in the wrinkled phase at b<b c (/j,,a) just like 
a at the continuous transition seen in Fig. 5(a). It should be noted that the 
scaling of a was not seen at the discontinuous transition point, where physical 
quantities might be ill-defined. 




Fig. 6. Snapshots of surfaces at (a) 6=1.86, a = — 5.5, L = 1.5Lq and at (b) 6 = 1.86, 
a = — 5.5, L = 3Lq. 

Snapshots of surfaces are shown in Figs. 6(a) and 6(b), which are obtained at 
b— 1.86, a = —5.5. The size is iV~ 1500 in each surface. The surface in Fig. 
6(a) looks short and swollen, whereas the surface in Fig. 6(b) long and slight, 
as expected. 
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5 Summary and conclusion 



To summarize the results, we have studied the scaling properties of string ten- 
sion of elastic membranes by grand canonical MC. The model contains a mea- 
sure term —a J2i logo"i in the Hamiltonian. The parameter a was assumed as a 
continuous one and fixed to a— — 5.5. The chemical potential \x were chosen so 
that iV~500, iV^lOOO, and iV~i500. Spherical surfaces were stretched, and 
then two boundary vertices were separated by length L— l.5L (N), 3L (N), 
where L (N)(cx y/~N) is the diameter of the sphere for the starting configu- 
ration of MC. The model undergoes the first-order phase transition on the 
surfaces of L = 3Lq(N) at a = —5.5. It was found that a becomes finite in the 
smooth phase at b>b c (/i, a) and vanishes in the wrinkled phase at b<b c (/i, a). 
When the length L is reduced to L — 1.5L (N), the phase transition changes 
to the second-order one, and a vanishes at the vicinity of the critical point of 
the continuous transition in the limit iV — ► oo as expected. 

We have explored the scaling of a at a — —5.5,0.0,5.5 with L = 1.5L (N), 
2L (N),3L (N), including the case presented in this Letter. Non- vanishing 
a was found only in the case a = —5.5 with L = 3Lo(N). More detailed 
information on the simulation data will be presented elsewhere. 

This work is supported in part by the Grant-in-Aid for Scientific Research 
15560160. 

References 

[1] F. David, in: D. Nelson, T. Piran, S. Weinberg (Eds.), Two dimensional 
quantum gravity and random surfaces, Vol.8, (World Scientific, Singapore, 
1989), p.87. 

[2] K. Wiese, in: C.Domb, J.Lebowitz (Eds.), Phase Transitions and Critical 
Phenomena, Vol. 19, (Academic Press, London, 2000), p. 253. 

[3] M. Bowick, A. Travesset, Phys. Rep. 344, 255 (2001). 

[4] J.F. Wheater, J. Phys. A Math. Gen. 27, 3323 (1994). 

[5] D.Sornette and N.Ostrowsky, in Micelles, Membranes, Microemulsions, and 
Monolayers, ed. W.M.Gelbart, A.Ben-Shaul and D. Roux, (Springer- Verlag, 
New York, 1994), p.251. 

[6] W. Helfrich, Z. Naturforsch, 28c, 693 (1973). 

[7] A.M.Polyakov, Nucl. Phys. B 268, 406 (1986). 



9 



[8] H.Kleinert, Phys. Lett. B 174, 335 (1986). 

[9] L. Peliti, S. Leibler, Phys. Rev. Lett. 54, (1985) 1690. 

[10] F.David and E.Guitter, Europhys. Lett. 5(8), 709 (1988). 

[11] F. David, Europhys. Lett. 52 (8), 577 (1986). 

[12] M.E.S.Borelli, H.Kleinert and A.M.J. Schakel, Phys. Lett. A 267, 201 (2000). 

[13] M.E.S. Borelli and H.Kleinert, Europhys. Lett. 53, 551 (2001). 

[14] H.Kleinert, Euro. Phys. J. B 9, 651 (1999). 

[15] Y.Kantor and D.R.Nelson, Phys. Rev. A 36, 4020 (1987). 

[16] Y.Kantor, in Statistical Mechanics of Membranes and Surfaces, Vol.5, ed. 
D.Nelson, T.Piran and S.Weinberg, (World Scientific, Singapore, 1989), p. 115. 

[17] G.Gompper and D.M.Kroll, Phys. Rev. E 51, 514 (1995). 

[18] J.F.Wheater, Nucl. Phys. B 458, 671 (1996). 

[19] M.Bowick, S.Catterall, M.Falcioni, G.Thorleifsson and K.Anagnostopoulos, J. 
Phys. I France 6, 1321 (1996); Nucl.Phys.Proc.Suppl. 47, 838 (1996); Nucl. 
Phys. Proc. Suppl. 53, 746 (1997). 

[20] H.Koibuchi and M.Yamada, Int. J. Mod. Phys. C 11(3), 441 (2000); 11(8), 
1509 (2000). 

[21] M.Bowick, P.Coddington, L.Han, G.Harris and E.Marinari, Nucl. Phys. Proc. 
Suppl. 30, 795 (1993); Nucl. Phys. B 394, 791 (1993). 

[22] K.Anagnostopoulos, M.Bowick, P.Gottington, M.Falcioni, L.Han, G.Harris 
and E.Marinari, Phys. Lett. B 317, 102 (1993). 

[23] S.M.Catterall, Phys. Lett. B 220, 253 (1989). 

[24] H.Koibuchi, Phys. Lett. A 300, 586 (2002). 

[25] H. Koibuchi, N. Kusano, A. Nidaira, K. Suzuki, M.Yamada, Phys. Lett. A 
319, 44 (2003). 

[26] H. Koibuchi, N. Kusano, A. Nidaira, K. Suzuki, M.Yamada, Phys. Rev. E 69, 
in press. 

[27] J.Ambjorn, B. Durhuus, and T. Jonsson, Phys. Rev. Lett. 58, 2619 (1987). 

[28] J.Ambjorn, A.Irback, J.Jurkiewicz and B.Petersson, Nucl. Phys. B 393, 571 
(1993). 

[29] F. David, Nucl. Phys. B 257 [FS14], 543 (1985). 

[30] D.V. Boulatov, V.A. Kazakov, I.K. Kostov and A. A. Migdal, Nucl. Phys. B 
275 [FS17], 641 (1986). 

[31] J.Ambjorn, B.Durhuus and J.Frohlich, Nucl. Phys. B 257, 433 (1985). 

[32] K.Fujikawa and M.Ninomiya, Nucl. Phys. B 391, 675 (1993). 



10 



